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Abstract 



In this dissertation, extensive researches on the thermal and electrical properties 
of millimeter-long aligned multiwall carbon nanotubes (MWNTs) prepared by thermal de- 
composition of hydrocarbons have been presented. 

The thesis consists of six chapters. Chapter 1 is an introduction, in which the 
discovery, structure and physical properties (including theoretical and experimental) are 
outlined. In Chapter 2, the preparation techniques and the characterizations of the carbon 
nanotubes are described. 

In Chapter 3, by using a novel self-heating 3w method, the specific heat, diffusivity 
coefficient and thermal conductivity of MWNTs are for the first time reported. MWNTs of 
a few tens nm diameter is found to demonstrate a strikingly linear temperature-dependent 
specific heat over the entire temperature range measured (10-300 K). The results indi- 
cate that inter-wall coupling in MWNTs is rather weak compared with in its parent form, 
graphite, so that one can treat a MWNT as a few decoupled two-dimensional (2D) single 
wall tubules. The thermal conductivity k of MWNTs has a roughly linear temperature 
dependence above ~120 K, below ~120 K, k follows almost a quadratic T-dependence. The 
room-temperature amplitude of n is found to be low, indicating the existence of substantial 
amount of defects in the MWNTs prepared by chemical-vapor-deposition method. 

In Chapter 4, the first true four-probe measurement of tunneling spectroscopy for 
single tunnel junctions formed between MWNTs and a normal metal is presented. The 
intrinsic Coulomb interactions in the MWNTs give rise to a strong zero-bias suppression 
of tunneling density of states (TDOS) that can be fitted numerically with the environmen- 
tal quantum-fluctuation (EQF) theory. At low temperatures, an asymmetric conductance 
anomaly near zero bias is observed, which is interpreted as Fano resonance in the strong 
tunneling regime. 

In Chapter 5, the thermoelectric power (TEP) and longitudinal magnetoresistance 
(MR) of MWNTs are measured. MWNTs have a moderately large positive TEP. At low 
temperatures, TEP has a metallic-like linear temperature dependence. The results give 
strong evidence that the electron-hole symmetry in metallic nanotubes is broken. The 
longitudinal MR of MWNTs is measured at 20 mK. Periodic oscillations are observed when 
a longitudinal magnetic field is applied. If only the outermost graphene wall of MWNTs 
contributes to electrical conductance, the period of oscillation agrees well with the period 
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h/2e of Altshuler-Aronov-Spivak (AAS) effect. The present results clearly indicate that 
quantum-interference dominate the resistance of MWNTs at low temperatures. 
Chapter 6 presents the main conclusions of this thesis. 



Keywords: graphite, carbon nanotubes, specific heat, thermal conductivity, 3a; method, 
tunneling spectroscopy, Coulomb interactions, Fano resonance, thermoelectric power, mag- 
netoresistance, AAS effect, quantum interference 
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Chapter 1 



Introduction 



1.1 Discovery of carbon nanotubes 

Carbon is really a special element in nature, whose chemical versatility makes 
it the central agent in most organic and inorganic chemistry as well as the fundamental 
constituent of all biological systems. 

Until recently, pure carbon was only known to exist in two forms: diamond and 
graphite. In 1985, Harold Kroto, Robert Curl and Richard E. Smally discovered a new 
form of carbon-the fullerenes, which are molecules of pure carbon atoms bonded together 
forming geometrically regular structures [1]. The best-known of these is Cqq, which has 
precisely the same geometry as a soccer ball, with one carbon atom at each point where 
the seams intersect. Because its structure is similar to that of the geodesic dome developed 
by the American architect, Buckminster Fuller, so this new molecule was named buckmin- 
sterfullerene, or buckyball. The discovery of fullerenes has added an important milestone 
in carbon research. 

Most importantly, research has proven the possibilities for producing a whole new 
range of carbon structures of various sizes, shapes and dimensionalities. According to Euller 
theorem, infinite kinds of closed cage structure can be made by add a minimum number 
of 12 pentagons to carbon honeycomb lattices. Soon later people discovered various such 
structures, such as CVo, C76, Cs4, Cqq, etc. Naturally, the possibility of existance of tube- 
sized fullerene molecules, i.e., carbon nanotubes, has also been proposed. 

In 1991, S. Iijima firstly observed the existence of carbon nanotube while studying 
the surface of carbon electrodes used in an electric arc-discharge appartus which had been 
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used to make fullerenes [2] . He found needle- like carbon whiskers with diameter of 4-30 nm 
and length of 1 /im. These carbon whiskers are hollow tubes made of concentric clinders 
of carbon honeycomb lattices. The distance between adjacent cylinders is similar to inter- 
layer distance of graphite, i.e., d ~ 0.34 nm. The discovery of carbon nanotubes stimulated 
tremendous interest in scientific community immediately. The subsequent success of the 
large-scale synthesis of nanotubes by T. W. Ebbesen and P. M. Ajayan of NEC opened the 
door for their widespread study. In 1993, Iijima [3] and D. S. Bethune [4] firstly observed 
single-wall carbon nanotubes (SWNTs) with diameter only ~ 1 nm. The relative simplicity 
of structure of SWNTs made them suitable for theoretical calculations. By 1996, A. Thess 
successly synthesized large quantities of high-quality SWNTs using pulsed laser evapora- 
tion technique [5]. Taking advantage of Scanning Probe Microscopy techniques obtainable 
nowadays, experimental explorations on the properties of individual SWNT made rapid 
progress, which has found extremely good agreement with theoretical predications. 

1.2 Structure 

Carbon nanotubes are made up of one or more wrapped-up seemless cylinders of 
carbon honeycomb lattices. The theoretically smallest nanotubes have a diameter equal to 
that of Cqo (d = 0.7 nm), which could be thought as two hemispheres of Cqo connected 
by a single wrapped-up graphite sheet. If a Cqq molecule is separated perpendicular to the 
quintic symmetry axis, the derived tube is called 'armchair' (Fig. |1.1| a); if it is separated 
perpendicular to the cubic symmetry axis, the derived tube is called 'zigzag' (Fig. |l.l] b). 
We can see that the orientation of carbon hexagonal rings is either parallel or vertical to the 
tube axis for these two kinds of tubes. Actually for most of the nanotubes, the orientation 
of hexagonal rings is neither parallel nor vertical to the tube axis, but has some fixed helical 
angel. Such nanotubes have chirality (Fig. |l.l] c). 



In 1992, Hamada et al. developed an indexing method to represent the structure 



for single shell of carbon nanotubes [7]. As shown in Fig. 1.2 a, cabon atoms are located at 
the corners of honeycomb lattice. The two primitive lattice vectors are ai,2 = ( a /2)(±l, \/3), 
where a = \/3d, with d = 0.1421 nm is the nearest-neighbour carbon separation. All the 
possible structures of carbon nanotubes can be represented by a pair of integers (n, m). A (n, 
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Figure 1.1: Schematic of the formation of carbon nanotubes by wrapping up graphene sheet. 
Different means of wrapping produce different tube structures: (a) armchair, (b) zigzag, (c) 
chiral. From Ref. [6] 
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m) nanotube is made by folding the two-dimensional hexagonal lattice shown and matching 
the lattice point (n, m) with the origin, with a wrapping vector C = na\ + ma,2- If we 
define the wrapping orientation of zigzag nanotube as the reference axis, the separation 
angle between the wrapping vector C and the reference axis is the helical angle 9. All 
nanotubes are uniquely defined by the helical angle 9 and the distance C which defines the 
tube size. Because of the symmetry of graphite lattice, 9 is chosen as < 9 < 30° to avoid 
duplicity. The three kinds of nanotubes have different helical angle 9. For armchair tubes, 
n = m, 9 = 30°. The zigzag tubes correspond to m = and 9 = 0. When n ^ m and not 
equal to zero, they are helical or chiral tubes with < 9 < 30°. 

For a (n, m) carbon nanotube, the diameter d and helical angle 9 are given by 

d= a c _ c (n 2 + nm + m 2 )^ (1-1) 

9 = arctan v3/(m + 2n) (1-2) 



1.3 Physical properties, theoretical 

Carbon nanotubes as a novel quasi one-dimensional material has stimulated great 
interest. From theoretical point of view, singlewall carbon nanotubes is ideal for study 
because of their relative simplicity. A straight-forward tight-binding calculation is based on 
the zone-folding picture, which deduces the band structure of carbon nanotubes from their 
parent form, graphene sheet, with periodic boundary conditions around circumference of 
nanotubes [7—11]. For a graphene sheet, each atom has one conduction electron in the 2p z 
state. The first Brillouin zone (BZ) is a hexagon whose sides are distant l/3a c _ c from its 
center. It is easily shown that this zone contains one electron per atom, so graphite is a 



half-filling semiconductor with zero gap, or a semimetal (shown in Fig. 1.3). Fermi surface 
reduces into six points at the corner of hexagonal BZ. Near Fermi level, the dispersion 
relation is nearly linear [12]. 

The periodic boundary conditions around circumference of nanotube make the 
crystal momentum transverse to the tube axis quantized. The allowed electronic states are 
restricted to those in the graphene BZ that satisfy [9] 



k ■ C = 2irm 



(1.3) 
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Figure 1.2: Schematic of the indexing method developed by Hamada et al. for single shells 
of carbon nanotubes. For example, a (6,3) nanotube is made by folding the two-dimensional 
hexagonal lattice shown and matching the lattice point (6,3) with the origin. The directions 
along which lattice points would make armchair and zigzag nanotubes are indicated. From 
Ref. [6] 
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Figure 1.3: (a) Energy dispersion of the 7r-band of graphene sheet near the K. The solid 
lines correspond to ID bands obtained by the zone-folding. Symbols are located at the 
positions where Van Hove singularities occur in these ID bands. The inset depicts a three 
dimensional view of graphene it/it* bands, (b) Comparison of the electron density of states 
obtained from scanning tunneling microscopy experiment (upper curve) and tight-binding 
calculations (lower curves). From Ref. [13] 



with m an integer (Fig. 1.4). In terms of the two-dimensional BZ of graphene, 



the allowed states will lie along parallel lines separated by a spacing of 2tt/\C\. Metallic 
behavior happens whenever the gapless Fermi points lie on the allowed transverse quantized 
wave vectors. For (n, n) armchair nanotubes, the gapless modes are present at k y = 0, 
resulting in metallic behavior independent of n. For zigzag or chiral tubes with (n, m), 
metallic behavior happens only if the vectors of Fermi points satisfy Eq. (1.3), leading to 
the condition n — m = 3q, where q is an integer. For the nonmetallic case, the smallest band 
gap for the nanotube will occur at the nearest allowed point to K, the Fermi wavevector, 
with amplitude given by tight-binding Hiikel model [9] 



E g = V ppn d/R (1.4) 

where V ppn is hopping matrix element of the 7r-bands, R is the nanotube radius(i? = \C\/2tt). 

This simplified model does not consider the symmetry breaking introduced by 
curvature effects. Due to curvature effects, the Fermi points will be shifted slightly along k x , 
but leaves the armchair tube gapless (Fig. |1.4| i). For chiral tubes which satisfie n — m = 3q, 
Fermi points are shifted away from the allowed quantized wave vectors slightly by order 
of 1/R 2 . Thus chiral tubes are quasimetallic with second-order small band gaps (Fig. 
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(a) 



(b) 




Figure 1.4: Fermi points in the Brillouin zone of carbon nanotubes. Dark circles and crosses 
indicate the locations of the gapless points without curvature-induced symmetry breaking, 
while gray symbols schematically indicate the shifted positions for curvature effects. Dashed 
lines cut the zone at a discrete set of allowed transverse momenta in (a) the armchair tube 
and (b) the zig-zag tube (here with N=3). From Ref. [15] 



i. no [14,15]. 

From above discussion, the armchair nanotube could be treated as truly ID metal. 



The low-energy conduction modes occur near the two gapless Fermi points ±K 



2tt (2 
- a \3' 



0). 



Excitation of other transversal modes costs the energy ~ 10 eV/n, and hence a ID condition 
arises. For (10, 10) tubes, the conditions for the low-energy regime are met even at room 
temperature. Up to energy scales ~ 1 eV, the dispersion relation around the Fermi points 
remains nearly linear (Fig. |1.5| a) and cause a finite and constant density of states (DOS) 



near Ep(Fig. 1.5 b). 



Based on studies of the electronic structures of SWNTs, theoretists also carried 
out some elementary studies on electronic strutures of MWNTs [17-19]. It was argued that 
a two-layer MWNT constitued by two concentric metallic (or quasimetallic) SWNTs is still 
metallic (or quasimetallic). An interesting result is, when one of the metallic nanotube is 
replaced by a semiconducting nanotube, a concentric metal-semiconductor device could be 
constitute [20]. But there also exists result that indicating changes of electronic properties 
due to inter-layer interactions when two SWNTs joined together. 

Another important research area is on the Mechanical properties of carbon nan- 
otubes. It is well-known that the in-plane C-C covalent bond in graphite is one of the 
strongest bonds in nature. The best estimate of the in-plane elastic modulus is 1.06 TPa 
and that of the tensile stiffness is 0.8 TPa. The mechanical properties of ideal carbon nan- 
otubes with negligible number of structural defects should approach or even surpass the 
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Figure 1.5: (a) Schematic of the electronic structure of (10,10) armchair nanotubes. (b) 
Corresponding finite and constant DOS near Ep. The 'gap' associated with DOS peaks at 
next one-dimensional subbands. From Ref. [16] 



reference values of graphite. Moreover, the rolled-up structure of carbon nanotubes would 
make them withstand repeated bending, buckling, twisting and compression at high rates 
with complete elasticity. 

The hollow structure of carbon nanotubes make it possible to fill them with other 
materials to form one-dimensional nano- wires. Theoretical studies point out that only- 
materials with small enough surface tensity could be filled into carbon nanotubes, including 
all the gases and liquids, most of the organic materials and some of the metals and their 
compounds [21-23]. 



1.4 Physical properties, experimental 

Although numerious theoretical calculations have predicted many novel physical 
and chemical properties for carbon nanotubes, the nanometer size and random orientations 
of nanotube samples make it extremely difficult to examine and certify these properties 
experimentally. However, taking advantage of the rapid progress of nano-fabrication and 
nano-manipulation, scientists are making fast progress on experimental studies and many 
valuable results has been obtained which agree quite well with theoretical predictions. 

Since nanotubes are ideal model systems for the investigation of low-dimensional 
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molecular conductors, measuring the electronic properties of individual nanotubes is always 
the focus of experimental studies. This is very challenging for two reasons. Both high- 
quality nanotube samples and new techniques for making electrical contacts to individual 
tubes are necessary. The first to be studied was individual multiwall nanotubes. By using 
STM lithography technique, Langer et al. reported the first measurement on individual 
MWNTs [24] . They found that transport properties of MWNTs were consistent with the 
quantum transport behaviors of weakly disordered low-dimensional conductors. Ebbesen 
et al. systematically measured four-probe conductance of individual MWNTs and observed 
both metallic and non-metallic behaviors [25]. Typical temperature dependence of the re- 
sistance of metallic tubes was similar to that of disordered semi-metallic graphite, with a 
moderate negative temperature coefficient. Another elegant method of electrical contact 
was used by Frank et al. [26]. A single MWNT attached to a STM tip was repeatly im- 
mersed and pulled out of a liquid metal like mercury. Astonishingly, an almost universal 
quantized conductance at room temperature was measured, providing evidence that trans- 
port in MWNTs is ballistic over distances of the order of > 1 //m. Recently, a pronounced 
Aharonov-Bohm resistance oscillation was observed by Bachtold et al. [27]. Their results 
clearly demonstrate that at low temperatures quantum interference phenomena dominate 
the magnetoresistance of MWNTs. 

Singlewall nanotubes have well defined structures and relatively less defects. The 
success of synthesis of high-quality SWNTs with uniformed structures greatly stimulated 
experimental studies. The first results on individual SWNTs were obtained by Tans et 
al. [28]. At low temperatures step-like current-voltage characteristics were observed that 
indicate single-electron transport with Coulomb blockade and resonant tunneling through 
single molecular orbitals. Quite remarkably, the nanotubes appear to behave as coherent 
quantum wires. The density of states of the molecule consist of well-separated discrete 
electron states corresponding to an one-dimensional 'particle-in-a-box' situation induced 
by micrometer tube length. The measurement on individual SWNT ropes by Bockrath et 
al. [29] obtained similar results. Soon later, two research groups at Delft Univ. and Harvard 
Univ. obtained the first atomic resolution STM images of SWNTs and observed the chiral 
winding of the hexagons along the tube [30,31]. The STM was also used to measure the 
tunneling spectroscopy of nanotubes. They found both metallic and semiconducting spectra. 
Their data provided the first experimental verification of the bandstructure predictions. The 
observed gap structures quantitatively agree with the calculations. 
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The explanations of all these experiments were based on single-electron picture, 
wherein the Landau Fermi liquid theorem is appropriate. However, the low-dimensional 
nature of carbon nanotubes would probably introduce pronounced electron-electron inter- 
actions, and a non- Fermi liquid behavior such as a Luttinger Liquid is predicted [32-35]. 
The first generation devices used in transport measurement of individual SWNTs have very 
high contact resistance, therefore the effect of electron-electron correlations was masked by 
charging effects. Tans et al. exhibited the first evidence of electron correlations by finding 
that electrons tunneling into a SWNT are spin polarized [36]. By improving the quality 
of electric contact, Bockrath et al. recently observed a power-law suppression of tunneling 
DOS on individual SWNT ropes [37], which could be well explained by tunneling into the 
bulk and ends of a clear Luttinger Liquid. 

The experimental characterizations of mechanical properties of carbon nanotubes 
is very difficult because of the short length (1-100 /zm) of the samples. Up to now, there 
are no direct tensile tests available to determine the axial strength of nanotubes. By mea- 
suring the amplitudes of intrinsic thermal vibrations of nanotubes during imaging inside a 
Transmission Electron Microscope (TEM), Treacy et al. estimated a Young's modulus as 
high as 1.8 TPa which is higher than the in-plane modulus for single-crystal graphite [38]. 
Recently, Pan et al. reported the first direct tensile test on MWNTs [39]. In their exper- 
iment, a tiny stress-strain puller was used to pull a MWNT bundle containing several ten 
thousand millimeter- long aligned nanotubes, and the average Young's modulus and tensile 
strength obtained were 0.45 ± 0.23 TPa and 1.72 ± 0.64 GPa, respectively. The relatively 
lower Young's modulus was attributed to the existance of large amounts of defects in CVD 
derived nanotubes, but it is still higher than that of Steel. Moreover, TEM observations 
have shown that carbon nanotubes could sustain extreme strain without showing signs of 
brittleness, plastic deformation or atomic rearrangements and bond rupture [40,41]. The 
flexibility is related to the ability of the carbon atoms to rehybridize. All these results 
indicate that nanotubes do have exceptional mechanical properties. 

1.5 Motivation and goal 

Despite the fast progress in the researches of carbon nanotubes, many important 
physical properties, such as thermal properties, have not been well experimentally char- 
acterized. Success of the synthesis of millimeter-long aligned multiwall carbon nanotubes 



Chapter 1: Introduction 



11 



opened the door of macroscopic manipulations on these nano-sized materials [42]. The 
continusity and good alignment of nanotubes enable us to perform measurements on well- 
defined samples. This dissertation aims at a systematic investigation on MWNTs, including 
their thermal properties, tunneling spectroscopy and electrical transport properties. 

On thermal properties, our goal is to explore the dimensionality of the phonon 
structure and its relationship with tube diameter. A potential rewarding experiment would 
be to measure the specific heat and thermal conductivity of a small sample simultaneously. 
Fortunately, a so-called 3a; method, which uses a narrow-band technique in detection, and 
therefore gives a relatively better signal-to-noise ratio, is especially suitable for such mea- 
surements. To do this, an explicit solution for the ID heat-conduction equation of a rod- 
or filament-like specimen is needed. 

On tunneling spectroscopy, remarkably results have been obtained on individual 
SWNTs, such as chirality related gap structure, one-electron tunneling, Luttinger liquid 
behavior, etc. However, no systematic measurement has been performed on MWNTs by 
now. MWNTs have diameters of several tens of nm. This is one order of magnitude higher 
than that of SWNTs and relates into an order of magnitude difference in energy scale. In 
zone-folding picture, MWNTs may still be treated as ID conductor, then an interesting 
and important problem is whether or not electron correlations exist in MWNTs. Taking 
advantage of the very long aligned MWNTs, we can use true four-probe configuration to 
measure the tunneling spectroscopy on a thin bundle sample containing only several tens 
of nanotubes. Important information will be obtained on the possible effects of electron 
correlations and their relations with tube diameter. 

Finally, electrical transport properties such as thermoelectric power, magnetore- 
sistance will be investigated. For thermoelectric power measurement, the continusity of 
the sample guarantees the reliability of the results. For magnetoresistance measurement, 
the larger transverse size of MWNTs compared with SWNTs make them more suitable to 
explore the orbital effects. A magnetic field of 3T is enough to induce a magnetic field flux 
of h/2e in a MWNT with outer diameter 30nm. It will be interesting to perform magneto- 
transport experiments of MWNTs in parallel and perpendicular field at low temperatures 
to explore quantum transport phenomena. 
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Chapter 2 



Sample growth and 
characterization 

2.1 Introduction 

For the first time, millimeter- long aligned multiwall carbon nanotubes were syn- 
thesized using thermal decomposition method by Dr. Z. W. Pan in Prof. S. S. Xie's research 
group [1]. Here a brief introduction of the sample growth and characterization is presented 
to let the readers have a first impression on the structure and morphology of MWNTs stud- 
ied in this dissertation. All the materials in this chapter were kindly provided by Prof. S. 
S. Xie. 

2.2 Process of synthesis 

The method used is similar to that reported by W. Z. Li et al. [2], but some 
improvement has been made in substrate preparation. Instead of using bulk mesoporous 
silica containing iron nanoparticles embedded in the pores as substrate, silica film on the 
surface of which iron/silica nanocomposite particles are evenly positioned was utilized. The 
iron/silica substrate was prepared by a sol- gel process using the technique described in Ref. 
[1]. Tetraethoxysilane (10 ml) was mixed with 1.5 M iron nitrate aqueous solution (15 ml) 
and ethanol (10 ml) by magnetic stirring for 20 min. A few drops of concentrated hydrogen 
fluoride (0.4 ml) was then added, and the mixture was stirred for another 20 min. The 
mixture was then dropped onto a quartz plate to form a film of thickness 30-50 //m. After 
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Figure 2.1: SEM image of the iron/silica substrate after reduction, showing iron/silica 
nanocomposite particles with size of 5-50 nm present evenly on every surface of the 
substrate. 



gelation of the mixture, the gel was dried overnight at 80° C to remove the excess water 
and other solvent, during which the gel cracked into small pieces of substrates of area 5-20 
mm 2 . 

The substrates were placed in a quartz boat and were then introduced into a 
chamber of a tube furnace. The substrate were calcined at 450°C for 10 h under vacuum 
(< 10~ 3 Torr) and then reduced at 500°C for 5 h in a flow of 9% hydrogen in nitrogen under 
180 Torr. At this stage, large quantities of nanoparticles with size of 5-50 nm formed evenly 



on all surfaces of the substrates (Fig. 2.1). Energy-dispersive X-ray (EDX) spectra taken 
from these nanoparticles showed the presence of iron, silicon and oxygen (68.7, 14.5 and 16.8 
wt%, respectively), which indicated that these nanoparticles were iron/silica nanocomposite 
particles. The iron/silica particles were believed to act as catalysts for nanotube growth. 
Subsequently, a flow of 9% acetylene in nitrogen was introduced into the chamber at a flow 
rate of 110 cm 3 /min, and carbon nanotubes were formed on the substrate by deposition 



of carbon atoms from decomposition of acetylene at 600° C under 180 Torr (Fig. 2.2). The 
growth time varied from 1 to 48 h. 

The as-grown nanotubes were examined by a scanning electron microscope (SEM; 
S-4200, Hitachi), and EDX spectra were recorded by a SiLi detector attached to the SEM. 
A transmission electron microscope (TEM; JEOL JEM-200cx at 200 keV) was used to 
characterize the fine structures of nanotubes. 
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Figure 2.2: Low-magnification SEM images of a sample after 48 hours of growth, (a) Only 
one nanotube array remains after stripping off other arrays from the other sides of the 
substrate, (b) Two nanotube arrays grown from the upper and lower sides of the substrate. 
From Ref. [1] 

2.3 Characterization 



SEM and TEM 

High-magnification SEM image (Fig. |2.3^ t) shows that carbon nanotubes grow out 
separately and perpendicularly from the substrate to form an array. The nanotubes within 
the array are of uniform external diameter (20-40 nm) and spacing (~ 100 nm) between 



tubes (Fig. 2.3 3). Most of the nanotubes in the array are highly aligned, although a few of 
them appear to be slightly tangled or curved. We note that no traces of polyhedral particles 
or other graphitic nanostructures are detected in both the bottom part and central part of 
the array, which indicates that the nanotubes prepared in this study have very high purity 
and thus have very high quality. 

The nanotube array can be easily stripped off from the substrate without destroy- 
ing the array's integrity, and the SEM images taken from the bottom end of the array 



confirm that the nanotubes are highly aligned and well separated (Fig. 2.3:, d). EDX spec- 



tra collected from the bottom end of the array demonstrate the presence of carbon alone, 
neither silicon nor iron could be detected. 

Low-resolution TEM observations of the bottom end of the nanotube array reveal 
that no seamless caps exist at the bottom ends of the nanotubes and the tubes are indeed 
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Figure 2.3: High-magnification SEM images of aligned carbon nanotube arrays, (a) Aligned 
nanotubes grow outwards perpendicularly from the surface of the substrate to form nan- 
otube arrays, (b) Higher magnification of the nanotube array showing aligned nanotubes 
have uniform external diameters (20-40 nm) and spacing (~ 100 nm) between tubes, (c, 
d) Bottom structure of the nanotube array showing well-separated, aligned and open nan- 
otubes. From Ref. [3] 
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Figure 2.4: (a) Low-resolution TEM image of the bottom end of the nanotube array show- 
ing aligned open carbon nanotubes. Encapsulated iron/silica nanocomposite particles are 
absent, (b) High-resolution TEM image of the bottom end of an open nanotube showing 
bowl-shape structure (see arrow) formed by the graphite sheets. From Ref. [3] 



opened (Fig. 2.4). The inner and outer diameters of the open nanotubes are 10-15 and 
20-40 nm, respectively. Encapsulated particles are not observed at the bottom ends of 
the nanotubes. High-resolution TEM observations confirm that the nanotubes are opened 
at the bottom ends. The nanotubes are well graphitized and typically consist of 10-30 
concentric graphite layers. 



Raman spectroscopy 

Raman spectroscopy is a very useful tool to characterize carbon materials because 
it gives information on the hybridization (sp n ) state, the size of the graphite crystallites and 
the degree of ordering of the material. Crystalline graphite, such as highly oriented pyrolitic 
graphite (HOPG), belongs to the space group Dg h with the vibrational modes having the 
irreducible representation T = 2E2 g + E u + Ayu + 2i?2g • Only the two E<i g modes are Raman 
active. These correspond to the movement of two neighbouring carbon atoms in a graphene 
sheet. HOPG is characterized by a single strong band at 1580 cm -1 which is one of these 
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£"2 9 modes. On the other hand, glassy carbon has its most intense band at approximately 
1350 cm -1 which is normally explained by the relaxation of the wave vector selection rule 
due to the finite size of the crystallites in the material. 

For raw MWNTs made by the arc-discharge method, the first-order Raman spectra 
has a weak and broad peak at 1346 cm -1 (D band) and a strong peak at 1574 cm -1 (G 
band) [4]. The G band is very similar to 1580 cm' 1 G band of HOPG but is slightly 
downshifted perhaps due to the curvature in the carbon network. The width of this peak 
is almost same to that of HOPG (approximately 23 cm -1 ), suggesting a highly crystalline 
graphitic structure. The 1346 cm -1 peak is a disorder-induced peak in graphitic structures 
and in nanotubes it could arise from the smaller crystallite sizes, defects in the curved 
graphene sheets or at the tube ends. Considering the Raman features, MWNTs must be 
seen as graphite-like microcrystals and are unlike fullerenes which show vibrational features 
of molecules. 

Z. W. Pan et al. studied micro-Raman spectroscopy of arrays of high-density 
aligned MWNTs made by thermal decomposition method and compared their results with 
that of HOPG and MWNTs made by arc-discharge method [5]. Here is a brief description 
of their results. 

Raman scattering experiments were performed at ambient conditions by a micro- 
Raman spectrometer (Reinshaw, England) in a backscattering geometry. The spectra were 
collected by using 514.5 nm excitation of argon laser at 0.2 mW and a spectral slitwidth 
< 2 cm -1 , with the laser incident perpendicular to the orientation of nanotube arrays. The 
diameter of focus area of laser beam is about 2 /xm. The time for collecting every spectrum 
is 5s x 3. 



Shown in Fig. 2.5 are Raman spectra of aligned MWNTs. Fig. 2.5 1 are first-order 
and second-order Raman spectra, Fig. |2.5b is third-order Raman spectrum. No smoothing 



was applied to the spectra. In Fig. 2.5 1, the 1570 cm -1 peak (G band) is the high-frequency 
Raman active E^g mode. It is downshifted by 11 cm -1 compared with G band (1581 cm -1 ) 
of HOPG, but is close to G band (1574 cm -1 ) of MWNTs made by arc-discharge method. 
Remarkablely, a strong peak can be seen at 1320 cm -1 (D band), and a weak peak at 1598 
cm -1 (D' band) appears. It is noticeable that the strong 1320 cm -1 peak resembles the 
D band of pyrolitic graphite (PG), but can not be observed in HOPG. In other carbon 
materials, the appearance of D and D' band is attributed to the crystallite-size effect or 
lattice twisting. The ratio of intensity of D band relative to G band reflects the size of 
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Figure 2.5: Micro-Raman spectra of aligned multiwall carbon nanotubes. (a) first-order 
and second-order Raman spectra, (b) Third-order Raman spectrum. From Ref. [5] 
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graphitized crystallites. In SEM and TEM studies, no carbon nano-particles was found. 
Therefore, the appearance of 1320 cm -1 peak could be attributed to defects in the curved 
graphene sheets or at the tube ends, turbostratic structure in graphene sheets, limited 
size of crystal domains in nanotubes and the coating of amorphous carbon at nanotube 
surfaces. Compared with MWNTs made by arc-discharge method, the intensity of D band 
is much higher. This mainly due to the much lower growth temperature of pyrolysis method 
(600° C). The results of Raman spectroscopy show that MWNTs made by pyrolysis method 
have relatively lower degrees of graphitization and larger amounts of defects. 

In the second-order Raman spectrum, four peaks at 2645 cm -1 , 2893 cm -1 , 2483 
cm -1 and 3219 cm -1 can bee seen. The 2645 cm -1 band is the second-order spectrum of 
D band (~ 2 x 1320 cm -1 ), the 2893 cm -1 band is identified as composition of D band and 
G band (1320 cm -1 + 1570 cm -1 ). 

In the third-order Raman spectrum, a weak 4220 cm -1 peak can be recognized. 
This peak is the third-order Raman spectrum of nanotubes and can be identified as compo- 
sition of D band and G band (2 x 1320 cm -1 + 1570 cm ) . It is downshifted compared with 
4320 cm -1 peak of HOPG and 4305 cm -1 peak of PG. No fourth-order Raman spectrum 
was observed in this experiment. 

From the above analysis, the Raman feature of MWNTs prepared has a strong 
resemblance to that of PG. This indicates that many nanometer-sized crystal domains exist 
in microscopic structure of MWNTs. The appearance of the third-order Raman spectrum 
shows that these crystal domains have high degrees of crystallization, which agrees with the 
conclusion from TEM studies. 

2.4 Summary 

Aligned and open-ended multiwall carbon nanotubes (MWNTs) were prepared by 
Dr. Z. W. Pan at a very high yield by pyrolysis of acetylene over film-like iron/silica sub- 
strates. These nanotubes grow outwards separately and perpendicularly from the surfaces 
of the substrates to form aligned nanotube arrays. The nanotubes within the arrays are 
of uniform outer diameters (20-40 nm), with a spacing of ~ 100 nm between the tubes. 
After continuous growth of 48 hours, the lengths of nanotubes can reach 2-3 mm, which is 
1-2 orders higher than previously reported values. TEM studies indicate a high degree of 
graphitization. The multiwall nanotubes are composed of 5-30 concentric graphene cylin- 
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ders with average inter-wall distance of 0.34 nm. The inner diameters of nanotubes are 
3-10 nm. Features of micro-Raman spectra show a strong resemblance to that of Pyrolytic 
Graphite. The appearance of a strong D band indicates the existance of nanometer-sized 
crystal domains in MWNTs. The appearance of the third-order Raman spectrum shows 
that MWNTs have high degrees of crystallization. 
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Chapter 3 



Thermal properties 

3.1 Introduction 

While the electronic properties of carbon nanotubes have been intensively studied 
and shown to be peculiar, such as the sensitive dependence of the electronic structure on 
the diameter and chirality of the tubule, the thermal properties of the material are not well 
characterized experimentally, though they are as important as the former both from the 
viewpoint of basic research and for possible application. Benedict et al. [1] predicted that the 
temperature dependence of specific heat would depend on the diameter of the tubule, and 
might even show a sign of dimensional cross-over as the temperature is varied. However, 
an experimental test of the thermal properties of the material is frustrated because the 
samples usually available are in a mat-like form, in which each individual tubule, or rope of 
tubes, is rather short and only loosely connected to each other, which prevents a reliable and 
equilibrium thermal measurement. The recent success in synthesis of millimeter- long highly 
aligned array of M WNTs [2] enables us to perform specific heat and thermal conductivity 
measurements on well-defined samples. 

3.2 Experimental 

The MWNTs were synthesized through chemical vapor deposition (CVD) [2] . High 
magnification scanning electron microscopy (SEM) analysis shows that the tubules grow 
out perpendicularly from the substrate and are evenly spaced at an averaged inter-tubule 
distance of ~ 100 nm, forming a highly aligned array. High resolution transmission electron 
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microscopy (HRTEM) study shows that most of the tubules are within a diameter range 
of 20-40 nm. The mean external diameter is ~30 nm. A tubule may contain 10-30 walls, 
depending on its external diameter. The continuity of the tubules, which is of crucial 
importance for the interpretation of the thermal conductivity data, is indicated by SEM 
studies, and by Young's modulus measurement where several tenth of TPa, a value higher 
than that of steel, was reached [3]. The final samples used in our measurements were 1-2 
mm long bundles of MWNTs stripped off from the bulk array. The apparent cross-section 
of the samples spreads from 10 -10 to 10 -8 m 2 . The filling factor of MWNTs in the bundles 
is ~1.5%, estimated from SEM observation. This estimation may bring some uncertainty 
in determining the absolute values of the thermal conductivity and specific heat of the 
MWNTs. However, their temperature dependence is not affected by this uncertainty. 

Measurements of thermal conductivity and specific heat on such tiny but well 
defined samples were made possible by using a self-heating 3a; method. The basic idea of 
the method can be traced back to 1910 [4], and was carefully tested later on. If both voltage 
contacts of the sample are ideally heat sunk to the substrate, but keeping the in-between 



part suspended (illustrated in Fig. 3.1), an ac current of the form Iosinut passing through 
the sample will create a temperature fluctuation at 2u), which will further cause 3oj voltage 
harmonics, V^, across the voltage contacts. V&j can be solved explicitly to an accuracy of 
~1/81 [5]: 

2IqLR (dR/dT) , fQ 

tancfio = 2il>7 (3-2) 

where 7 = 1? /ir 2 a 2 , a 2 = n/C p p m is the diffusivity coefficient, k, the thermal conductivity, 
C p , the specific heat, p m , the density, R, the resistance, L, the length (between voltage 
contacts), and S, the cross-section of the sample. By measuring the frequency dependencies 
of the amplitude and the phase shift cpo of V&j in proper u and Iq ranges (typical results 



are shown in Fig. 3.1), both k, a 2 , and hence C p of the sample can be determined. This 
method was verified to be reliable on pure Pt-wire samples. The measured C p of Pt agrees 
with the standard data within an accuracy of 5% from 10 K to 300 K. 
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f (Hz) 

Figure 3.1: The frequency dependence of the amplitude (main frame) and the phase shift 
4>q (upper right inset) of the 3w voltage measured on one of the samples at T = 50 K. Solid 
lines are the predictions of Eqs. 1 and 2. Lower left inset: schematic diagram of sample 
mounting. The samples were glued to sapphire substrate with silver paint in such a way 
that the inner part between two voltage contacts was suspended. The silver paint contacts 
serve both as electrodes and heat sinkers of the sample. The 3w voltage was measured 
by using a lock- in amplifier (SR850). The whole set-up was heat-shielded to the substrate 
temperature and was maintained in high vacuum, to preclude heat leakage through air. The 
total heating power at each measured temperature was adjusted to be small comparing to 
the value of sample's thermal conductance. The heat leakage through radiation is estimated 
to be less than 5%. 
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3.3 Results and discussion 



Shown in Fig. 3.2 are the main results. First we notice that, while both k and a 2 
are significantly non-linear, C p calculated from them follows a pretty linear T-dependence 
over the entire temperature range measured. We note that a T-linear phonon specific heat 
was predicted for some special systems, but seldom observed experimentally over such a 
wide temperature range. This behavior is dominated by the phonon contribution, since the 
electron contribution is negligible in the temperature range of this experiment, as expected 
from the fact that the electron density of states is either very low or gapped at the Fermi 
surface of carbon nanotubes [6-8]. 

The phonon contribution can generally be written as: 

c„=r-* B (^y^^ ( 3 .3) 

where p(uj) is the phonon spectrum determined by the phonon energy dispersion u>(k) 
of different modes, and by the occupational dimensionality of phonon excitations in k- 
space. Detailed explanation of the T-dependence of C p requires the information on the 
p(co) of the MWNTs. While the p(co) and u(k) for single-wall tubules of ~1 nm diameter 
have been calculated numerically and shown to have a linear phonon dispersion and ID 
behavior [10-12], such calculation is virtually impossible to be performed on tubules of a 
few tens nm diameters, containing over hundreds of hexagonal cells along the circumference. 
Nevertheless, one can still find immediately from Eq. 3.3 that a constant p{oo) would result 
in a linear T-dependence of C p when T <§C @d = ^w^jj /&B, the Debye temperature. 
This condition is believed to hold in our case, referring to the fact that the C-C bond 
configuration within the wall is pretty much the same as in graphite, hence its Qd should 
not be far from ~2400 K [9]. Therefore, we conclude that the linear T-dependence of C p 
represents a constant p(co) of the MWNTs for the phonon states excitable in the temperature 
range of this experiment. It is less likely that a complicated p(ui) could "happen" to result 
in a linear C p after the integration, in this case one would usually get a complicated T- 
dependence of C p . Strictly speaking, due to the rolled- up nature of the material, there 
should exist some spiky- like structures over the flat background of MWNTs' p(uj). Such 
structure is most significantly seen in tubules with smallest diameter [11]. Nevertheless, 
thermal properties like C p are integration over a wide frequency range and hence would be 
rather insensitive to the local spiky structures. 
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T (K) T (K) 

Figure 3.2: (a) Temperature dependence of the thermal conductivity of three MWNT sam- 
ples. It is roughly linear above ~120 K, but becomes quadratic at temperatures below: 
k oc j ,1 - 98±003 (the solid line), (b) The thermal diffusivity coefficient of the samples, (c) 
The specific heat of the samples, calculated using the data in figures (a) and (b), showing a 
linear temperature dependence, which differs significantly from that of graphite [9] (see the 
inset). For graphite, C P (300K) ~ 720 Wm _1 K _1 . (d) Typical resistance R (main frame) 
and dR/dT (inset) of the samples, which are used by Eq. 3.1 to calculate the k and a 2 . 
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A constant p(uj) usually originates from linearly dispersed branches in a ID system, 
or from one or more quadratically dispersed branches in a 2D system. The dimensionality of 
our MWNTs can be determined through the following criteria [1]. Because of the rolled up 
structure, the allowable phonon states are quantized in fc-space, forming some discrete lines 
parallel to the tubule axis. A tubule is strictly ID-like if only those states on the central 
line are thermally excitable, otherwise if many lines are occupied with phonons then the 
tubule should be 2D-like. A simple estimation shows that even at the lowest temperature 
of this experiment (~10 K), and even for the inner-most wall of our MWNTs (which has 
the largest inter-line distance), the thermally excitable phonon states would cover several 
to several tens of lines, depending on whether a linear or a quadratic dispersion is assumed. 
Therefore, our MWNT at higher temperatures (i.e., >10 K) can approximately be regarded 
as a 2D system with a constant p(oj). This approximation sounds better if considering the 
fact that the outer walls, which have larger weight on the total C p than the inner walls, 
undergo dimensional cross-over at even lower temperatures. 

Such a phonon structure is similar to, but still different from that of graphite 
whose p(u) is only roughly flat in a moderate frequency range [13]. The similarity is easy to 
understand, for that graphite is the parent material of MWNTs. The central issue now is 
how to understand the difference in p(oj) induced by the rolling up of graphene sheets. For 
flat sheet, as in graphite, its hexagonal symmetry leads to a quadratically dispersed acoustic 
phonon branch, corresponding to the out-of-plane vibration mode of the sheet. This branch 



contributes to a constant term in p(oo) [13-16] (Fig. 3.3). However, there are at least two 
reasons that prevent graphite from strictly demonstrating a constant p(u) and a linear C p . 
First, inter-layer coupling could not be neglected at low frequencies, which gives p(uo) ~ w 
and C p ~ T 2 (in the temperature range of interest here we do not consider the case of very 
low frequencies). Second, the other two linearly dispersed in-plane branches, which gives 
p(io) ~ ui and C v ~ T 2 , also contribute to the total specific heat [14,15]. Therefore, no 
strictly T-linear C p as well as strictly flat p{oS) exist [9]. 

From the above discussion, the fact that MWNTs' linear C P (T) extrapolates to the 
(0, 0) point seems to indicate that the out-of-plane acoustic mode (as in a graphene sheet) 
dominates the thermal properties of this 2D-like material, and that the two in-plane acoustic 
modes, presumably linearly dispersed, give a negligible contribution. Here an assumption 
has spontaneously been made: for MWNT of a few tens nm diameter its acoustic phonon 
modes can still be classified as one LA and two TAs as in graphite. This assumption should 
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Figure 3.3: Graphite phonon dispersion relations and the corresponding phonon density of 
states. From Ref. [17] 



be true, for that at most temperatures of this experiment the wavelengths of the majority of 
phonons are much shorter than tubules' diameter, i.e., many lines in /c-space are occupied 
with phonons. The striking linearity of C p further indicates that the coupling between the 
walls in MWNTs is much weaker than that in graphite, so that one can treat a MWNT as 
a few decoupled single-wall tubules, as far as their vibrational properties are concerned. 

Two facts may be responsible for the weak inter- wall coupling: the larger inter-wall 
distance in MWNT than the inter-layer distance in graphite, and the turbostratic stacking 
of adjacent walls which is unavoidable in the rolled-up structures. Although there are a 
number of experiments indicating that the inter-wall distance in MWNT is larger than that 
in graphite [18,19], recent HRTEM study [20] shows that the inter- wall distance decreases as 
tubule's diameter increases, at diameters >10 nm it saturates to ~0.344 nm, a characteristic 
inter-layer distance in turbostratic stacking. Therefore, we believe that the weak inter-wall 
coupling in MWNT is rather caused by the turbostratic stacking of adjacent walls. 

The thermal conductivity of our MWNTs is rather low compared to what is gen- 
erally expected [21]. Since the apparent Wiedemann-Franz ratio R(kS/L)/T is about 
two orders of magnitude larger than the free electron Lorenz number, the measured n 
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essentially reflects the phonons contribution, similar to the cases in graphite and carbon 
fibers. Whereas k of graphite shows a marked maximum around 100 K, n of our MWNTs 
monotonously decreases with lowering temperature. If roughly expressing k as k ~ T n , 
then n is a little less then unity above ~120 K, resulting in a faint hump of diffusivity 
a 2 between 120-300 K. Below ~120 K, k follows almost a quadratic T-dependence, i.e., 
n ~ 1.98, leading to a quick decrease of a 2 , or to say, a nearly linear decrease of apparent 
phonon mean-free-path with temperature. This behavior indicates that drastic change in 
phonon scattering mechanism occurs at ~120 K. Comparing with graphite, although the n 
of the latter also roughly takes a quadratic temperature dependence at low temperatures, 
the mean-free-path of phonons in graphite is rather confined by boundary scattering, be- 
cause the C p of graphite is quadratically temperature-dependent at low temperatures (i.e., 
below 70 K [22]). An energy-dependent mean- free-path in MWNTs contradicts to the usual 
boundary scattering mechanism seen in graphite and other solids at low temperatures. One 
difference between nanotubes relative to graphite sheets is the existence of the radial breath- 
ing mode in the former [23]. Perhaps the thermal excitation of such phonons is responsible 
for the change in scattering mechanism around 120 K. Further investigation is needed to 
clarify the issue. 

The MWNTs grown by CVD method at temperatures as low as ~600 K are far 
from perfect, as indicated by the fact that their thermal conductivity and electrical conduc- 
tivity are about two orders of magnitude lower than those of perfect crystalline graphite at 
room-temperature. The average electrical resistivity of MWNTs was estimated by assum- 
ing that all walls in a MWNT conduct electric current, the characteristic value is about 
1 x 10~ 4 ilm at 300 K. This value, together with the weak negative temperature-dependent 
coefficient, are similar to those of disordered semimetallic graphite, and are comparable to 



previously reported data on MWNT bundles [25-27]. Figure 3.4 shows that the T-dependent 
behavior of MWNTs resembles those of not well graphitized materials such as the as-grown 
carbon fibers and glassy carbons [24] , whose in-plane domain size of the hexagonal ordering 
is only a few nm. It has been shown in those materials that the mean-free-path of the 
phonons deduced from the k data is well correlated with the domain size [24,28]. We there- 
fore speculate that the domain size of our CVD-derived MWNTs is the same as the phonon 
mean- free-path I, which is a few nm if estimated using the classical relation k ~ C p vl and 
assuming a characteristic sound velocity v ~ 10 4 m/s. This picture is further supported 
by micro- Raman investigation wherein the 1320 cm" 1 band, which is prohibited in perfect 
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10 100 
T (K) 

Figure 3.4: A log-log plot of the thermal conductivity of our samples (circles) comparing 
with that of vapor-grown carbon fibers [24] (triangles: as-grown, squares: heat treated at 
3000 °C). Inset: First-order Raman spectrum of our MWNTs. 
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graphite, becomes substantially high in our MWNTs, resembling the case of as-grown CVD 
carbon fibers [24]. Nevertheless, the gross phonon spectrum is unlikely to be significantly 
altered upon the reduction of the domain size to a few nm [24]. We also note that for 
disordered systems, a two-level mechanism could contribute to a linear term of specific heat 
additive to the ordinary term and even surpass the latter at very low temperatures [29]. 
Obviously, such a mechanism, if exists, cannot account for a total linear C p of the MWNT 
over the entire temperature range measured. Therefore, our previous discussions on the C p 
and the phonon structure of the MWNTs should still be valid. 

3.4 Summary 

The specific heat and thermal conductivity of millimeter-long aligned carbon mul- 
tiwall nanotubes (MWNTs) have been measured. As a rolled-up version of graphene sheets, 
MWNT of a few tens nm diameter is found to demonstrate a strikingly linear temperature- 
dependent specific heat over the entire temperature range measured (10—300 K). The results 
indicate that inter-wall coupling in MWNT is rather weak compared with in its parent form 
— graphite, so that one can treat a MWNT as a few decoupled two-dimensional (2D) single 
wall tubules. The thermal conductivity is found to be low, indicating the existence of sub- 
stantial amount of defects in the MWNTs prepared by chemical-vapor-deposition method. 



3.5 Appendix: 3£1 method 

We present here a novel and simple way to measure the specific heat and thermal 
conductivity of rod- or filament-like specimens. This method uses four-probe configuration 
same as for measuring electrical resistance, except for that the voltage contacts serves not 
only as electrical contacts but also as thermal contacts, to heat-sink the specimen at these 
two points to the substrate temperature. The specimen, which needs to be electrically 
conductive and with a temperature-dependent resistance, serves both as heater and tem- 
perature sensor to tell the thermal response of the material. With this method we have 
successfully reproduced the specific heat and thermal conductivity of Pt-wire specimens, 
and measured the specific heat and thermal conductivity of tiny carbon nanotube bundles, 
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some of which are only ~10 -9 g in mass. 

Specific heat and thermal conductivity are two of the most fundamental properties 
of the condensed matter. Over the past centuries, many methods have been developed for 
measuring these quantities. One noticeable category of the methods is the so called "3a/' 
method, which uses a narrow-band technique in detection, and therefore gives a relatively 
better signal-to-noise ratio, especially for the measurement on small specimens. In this 
method, either the specimen itself serves as a heater and at the mean time a temperature 
detector, if the specimen is electrically conductive and with a temperature-dependent elec- 
trical conductivity, or, for electrically non-conductive specimens, a metal strip is artificially 
deposited on the surface of the specimen to serve as the heater and detector. An ac current 
of the form Iosmcot passing through the specimen or the metal strip creates a temperature 
fluctuation at an angular frequency of 2w, and consequently, a resistance fluctuation at 
2cj, which further generates a voltage fluctuation at 3oj. Corbino [4] is probably the first 
person noticed that the temperature fluctuation of an ac heated wire gives useful informa- 
tion about the thermal properties of the constituent material. Systematic investigations to 
the 3cj method were carried out mainly during the 1960's [30-32], and in the recent ten 
years [33-38] which made this method a beautiful and practicable method. Nevertheless, 
we notice that in the previous works the heat-conduction equation were solved under the 
approximations either only for high-frequency limit [30,31,38], or only for low-frequency 
limit [33,35,36]. After these approximations one accordingly either loses the information of 
the thermal conductivity, or the specific heat of the specimen. 

Here we present an explicit solution for the ID heat-conduction equation. With 
this solution and together with the use of digital lock- in amplifier available nowadays, we 
are able to simultaneously obtain the specific heat and the thermal conductivity of a rod- or 
filament-like specimen through measuring the 3w voltage using a simple four-probe config- 
uration. We have tested this method on Pt wire specimens. Correct values of specific heat, 
thermal conductivity, and Wiedemenn-Franz ratio, were obtained. With this method, we 
are able to measure the thermal properties of very tiny specimens such as carbon nanotubes, 
sometimes only 10 -9 g in mass. 



In Section 3.5.1 we will present an explicit solution of the heat-conduction equation. 



In Section 3.5.2 we will give the experimental test of the method on Pt wires (results on 
carbon nanotubes are already presented previously). A summary of the method will be in 
the last section. 
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Figure 3.5: Illustration of the four-probe configuration for measuring the specific heat and 
thermal conductivity of a rod- or filament-like specimen. 



3.5.1 Solution of the ID heat-conduction equation 
in the presence of an AC heating current 

We consider a uniform rod- or filament-like specimen in a four-probe measurement 



configuration as shown in Fig. 3.5 The two outside contacts are used to feed an ac current of 
the form iQsinujt into the specimen, and the two inside ones are used to measure the voltage 
signals, just as in a standard electrical resistance measurement. The main differences from 
a pure electrical resistance measurement are (i) the part of the specimen in-between the two 
voltage contacts has to be suspended, to allow temperature variation along the specimen, 
(ii) All the contacts should be highly thermal-conductive, to heat-sink the specimen at 
these points to sapphire substrate, (hi) The specimen has to be maintained in high vacuum 
and the whole setup be heat-shielded to the substrate temperature, to preclude heat loss 
through air conduction and radiation. In such a configuration, the heat generation and 
diffusion along the specimen can be described by the following partial differential equation 
and the initial and boundary conditions: 



d 

pCp-Q-Tixit) 



d 2 



lQs'm 2 u>t 



[R + & (T(x,t)-T )] 



(3.4) 
(3.5) 



LS 

T(0,t) = T 
T(L,t) = T 
{ T(x,0)=T 

where C p , k, R and p are the specific heat, thermal conductivity, electric resistance and 
mass density of the specimen, respectively, R' = j& at substrate temperature To, L, the 
length (between voltage contacts), and S, the cross section of the specimen. 

Setting T(x, t)—To to T(x, t) will not affect the solution but simplifies the equation 
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and the initial and boundary conditions: 



d d 2 

—T(x, t) — a 2 — —^T(x, t) — c sin 2 ujt T(x, t) = b sin 2 wi 
at aar 



T(0,t) = 
T(L,t) = 
I T(x,0) = 



(3.6) 
(3.7) 



where a 2 = is the thermal diffusivity, b 



I%R 



l£_R' 



P Cp 10 "^^^ ^XLtOlVlUJ, V — pCpLS , ^ — pCp LS- 

Using the impulse theorem, T(x, t) can be represented as the integral of the spec- 



imen's responses to the instant "force", the right-side term of (3.6), at each time interval: 

(3.8) 



where z(x, t;r) satisfies: 



T(x,t) = / z(x,t;T)dT 



dz 2 d 2 z 2 

c sin cot z = 



dt dx 2 
z(Q,t) = 
z(L,t) = 

z(x, t + 0) = b sin 2 WT 
Expressing z(x, t; r) with the Fourier series of the form: 



(3.9) 
(3.10) 



z(x, t; t) = y~] U n (t; r)sin 



nirx 



n=l 



L 



(3.11) 



and subtracting (3.11) into (3.9), we have 

^2 



£ 

n=l 



+ 



n 



7 



c sin cut U n 



mrx 



sm- 







(3.12) 



where 7 = (L/na) 2 . 

The term c sin 2 wi can be omitted if n 2 /^ 3> c, or equivalently 



< 1 



(3.13) 



n 2 ir 2 KS 

This condition is usually well satisfied. For example, in a typical measurement Io=10 mA, 



R'=0.1 Q/K, L=l mm, S=10~ 2 mm 2 , k=100 W/m K, the left side of (3.13) is ~ 10" 



After omitting the c sin 2 cjt term, the solution of the ordinary differential equation (3.12) is 
then: 



U n (t;r) = C n (r)e -> 



(t-r) 



(3.14) 
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where C n (r) can be determined using the initial condition in (3.10), together with the 

(3.15) 



relation Y°°_, ^t^sin^ = 1 for < x < L: 

/— m— 1 nn L 



nir 



Using (3.14) and (3.15), (3.11) becomes: 



z(x,t;r) = ^2 



sm- 



nnx 2b[l - (-1) 



n=l 



L 



■ 2 - — (t-r) 
8111 ujt e 1 y ' 



mr 



(3.16) 



Subtract (3.11), (3.16) into (3.8) and reset T(x,t) to T(x,t) — Tq, we have: 



T(x,t)-T = J^sin 



71=1 



nirx 7&[1 
~L 



n 3 7r 



cos2a;t + ^ sin2a;t 



1 + 



/ 2ary 
V n 



(3.17) 



where we have already dropped off the exponential term which should fastly decay out with 
time. 

The total resistance of the specimen at temperature To is: 

r-L 



R = y f {R + R'[T(x,t) -T ]}dx 
L Jo 



(3.18) 



Using (3.17) and the relation svaV^dx = [1 — (—1)™]^, resistance fluctuation 



with respect to R at To can be expressed as: 

» 7& g[l_(.l)»]2 

R fl = l^ ^2 

n=l 



cos2cjt + ^ sin2wt 
1 _0f 

i + (^) 2 



(3.18) 



Obviously, the n = 2 term automatically vanishes. If only accounting the n = 1 
term, which introduces a relative error of the order ~ 3 -4 , we have: 



Rfi = -2rtR' 



7T- 



1 



sin(2u;i + 



1 + (2wyY 



where 4>o satisfies: 



tan^o = (2^7) 1 



(3.20) 



(3.21) 



The according voltage fluctuation in the presence of the ac current is Rfilos'mut. 
After dropping off the dc term and re-define a phase constant cf) = \ — (f>o, ac voltage has 
the form: 



2IILRR' 



tt a kS^1 + (2w7) 2 



sin(3a;i — <fi) 



(3.22) 
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tan</> = 2^7 (3.23) 
If using rms values of voltage, V 3uJ , and rms values of current, /, as what lock-in 



gives, (3.22) becomes: 



AI 3 LRR' 



V 3uj = . (3.24) 

l^KSy/l + (2 W 7) 2 V ' 

The following alternative form makes it more clearly how the 3uj signal depends 
on the dimensions of the specimen: 

3U1 + (2w 7 ) 2 \ 

where p e is the electrical resistivity of the specimen, p' e 



(3.25) 



dT ■ 



If accordingly only taking the n = 1 term in (3.17), the temperature fluctuation 
along the specimen is: 



276 nirx 
Tft(x,t) = sin— 

7T L/ 



1 



sin(2u;t + (f>o] 
V / l + (2^7) 2 



(3.26) 



3.6 



where T9 l = ir/2jb. Obviously, temperature fluctuation 



(shown in Fig. 3.6 (a)), whereas the temperature along 



The result is illustrated in Fig. 
reaches the maximum as ojj - 
the specimen saturates to a time- independent distribution as W7 — > 00 (Fig. |3.6| (c)). 

At the high-frequency limit wy 3> 1, i.e., when the thermal wavelength A -C L, 



where A is defined as A = yf|j, (3.24) becomes: 

2I^RR' 



V 3ui 



(ujj > 1) 



(3.27) 



■K 2 uipC p LS 

which is the same as Holland's result [30] except for a slight difference in the numerical 
coefficients (i.e., ^ here vs. | there). At this limit, however, one losses the information 
of thermal conductivity of the specimen. Oppositely, if taking the low-frequency limit, 
wy C 1, or A > I, one will only get the thermal conductivity of the specimen, but 
loss the information of specific heat, similar to what happened in Gahill's treatment for a 
two-dimensional heat diffusion problem [36]. At the low- frequency limit, V 3u} is frequency- 
independent: 

V** = ^-i^ (*7 « 1) (3.28) 

In above derivation we have neglected the radial heat loss through radiation. Such 
heat loss per unit length from a cylindrical rod of diameter D can be expressed as: 



W(x, t) = iro-D [T 4 (x, t) - T 4 ] rj AttoDTq [T(x, t) - T ] 



(3.29) 
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Figure 3.6: Illustration of the temperature variation along the specimen under a driving ac 
current Iosinut. The range of the variation is marked as shadowed area, which reaches the 
maximum at the limit ujj — > and shrinks to a line as ujj — > oo. 
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Equation (3.6) and (3.7) can be re- written as: 
8 



Ot 



T(x,t) 



dx 2 



T(x, t) (g — c sin uit) T(x, t) = b sin ut 



T(0,t) -- 

T(L,t) 

T(x,0) 





= 
= 



(3.30) 



(3.31) 



where 9 = ^gf- 



(3.12) now becomes: 
dU n 



E 

71=1 



alt 



+ 



n 



7 



+ g — c sin ujt I U n 



sm- 



nirx 



(3.32) 



After replacing — + g with ^V, the final solution of (3.32) is the same as (3.22) 



and (3.23), except that now 7 in (3.22) and (3.23) should be replaced with 7'. In the 



approximation of only taking the n 
heat loss can be neglected if 



1 term, we have: 7' 



91 < 1 



!+ g7 - Obviously, radiation 



(3.33) 



16<tT 3 L 2 

For cylindrical rod, (3.33) becomes: — 2 ° n <C 1. Shortening the length or increasing the 



diameter of the rod will help to eliminate the effect of heat loss through radiation. 

Similar discussion can also be applied to account for the heat loss by remnant air. 
But for detailed calculation one needs to know the thermal resistance of the specimen at 
the surface. 



3.5.2 Experimental test 

The test of the above solution was carried out on two kinds of specimens: Pt 
wires and bundles of multiwall carbon nanotubes. The electrical resistance of the former 
has a positive temperature coefficient, whereas the latter a negative coefficient. Within the 
reasonable parameter ranges of frequency and current (we will explain it later) , we do find 
that the is proportional to I 3 / 'Wl + (2wj) 2 , and the phase angle tan</> is proportional to 
the angular frequency. For Pt specimen, the yielded specific heat and thermal conductivity, 
as well as the Wiedemenn-Franz ratio, do agree with the standard data over the entire 
temperature range measured. 

We use a digital lock-in amplifier SR850 to measure the 3a; voltage signal. The 



block diagram of measurement is shown in Fig. 3.7 We turned off all the filters of the 
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digital lock-in 
SR850 or 830 
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ac current 
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current 
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Figure 3.7: Block diagram of the measurement. We choose the digital lock-in amplifier 
SR850 or 830 to measure the 3w voltage. The reference-out of the lock-in is boosted to be 
an ac current by a simple circuit (lower panel), and is then applied to the specimen. The 
feed-back resistor R* should be those with very low temperature coefficient. 



lock-in during the measurement, and used the dc coupling input mode. Before measuring 
the 3a; signal, the phase of the lock-in amplifier was adjusted to zero according to the la; 
signal. We use a simple electronic circuit to booster the la; reference-out voltage into an ac 
current feeding to the specimen. The 3a; component of the current is 10 -4 or less comparing 
with the la; component, as checked by an HP89410A spectrum analyzer. Because the 3a; 
voltage signal is deeply buried in the lo> voltage signal of many order of magnitude higher 
in amplitude, high dynamic reservation of the lock-in is usually required. We keep the 
dynamic reservation unchanged relative to the total magnification of the lock-in amplifier 
during the entire measurement. 

There are two ways to perform the measurement. In the first way the substrate of 
the specimen is maintained to a few fixed temperatures. At each temperature the frequency 
dependence of the V^d is measured. In this way one can check whether the I 3 and the 



l/yl + {2uj~f) 2 dependencies of V&j , as well as the tan0 = 2ujj relation, are held. 

Because V&j oc I 3 , one will get much larger signal if using larger /. However, 
there are two reasons that prevent us from using inadequately large /. Firstly, excessive 
heat accumulation on the specimen will create significantly large temperature gradient at 



the silver paste contacts, which will violate the boundary condition in (3.5). Secondly, 
radiation heat loss will increase significantly as the temperature modulation is increased. 
In all the cases the expected relation oc I 3 will no longer be held. On the other hand, 
the relation will also be violated if using too small I, since in the latter case could be 
comparable to, or even smaller than, the spurious 3a; signal from the current and other 
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sources. In our measurement the total heating power at each measuring temperatures was 
maintained to be approximately equal to the value of the specimen's thermal conductivity, 
i.e., keeping the temperature modulation to be around 1 K along the specimen. 

In the second way of measurement, the temperature of the substrate is slowly 
ramped at a fixed rate, meanwhile the working frequency of the lock-in amplifier is switched 
between a few set values. The whole process, including the temperature ramping rate, the 
amplitude of the heating current, the frequency set-points, are all controlled by a personal 
computer. 

For Pt specimen, we chose a wire of diameter D = 20 /im and length L = 8 mm. 
We found that 7 of the specimen varied from 0.005 s" 1 at 10 K to ~ 0.2 s -1 at room 
temperature, so that the working frequencies were chosen to be between 1 to 80 Hz, in 
order to cover the most featured frequency range of the functional form 1 / y/l + {2uj^) 2 . 

Shown in Fig.|3.8|(a) is the current dependence of the at 22K, demonstrating an 



J 3 -law in a mediate current range. Figure 3.8 (b) and (c) show the frequency dependencies 



of the amplitude and the phase angle of the V&u signal at 22 K, respectively. It can be seen 
that the data (open circles) agree well with the predicted functional forms (solid lines). By 



fitting the data in Fig. 3.8 (b) to the functional form of (3.24), we get the information of 
thermal conductivity k (Fig. |3.8] (d) , open circles) and the time constant 7. Note that the 7 



3.8 



(e) and (f) as open circles. 



39] (the solid circles in Fig. 



3.8 



thus obtained is in fairly good agreement with what the phase angle reveals via Eq. (3.23). 
Thermal diffusivity a 2 and specific heat C p of the specimen can be obtained via relation 
7 = (L/-7ra) 2 and a 2 = The results are shown in Fig. 

C p thus obtained agrees well with the standard data for Pt 
(0)- 

We notice that the measured thermal conductivity of the Pt wire only show a small 
peak at low temperatures, far insignificant than what is expected for high purity Pt material. 
Such difference is possible since k depends largely on the purity, structural perfection, and 
even the dimensions of the specimen. We believed that the n data we obtained reflect the 
true thermal conductivity of our Pt specimen, since the Wiedemenn-Franz ratio deduced 
from the k and the electrical resistivity, or more directly, from the thermal conductance 
and the electrical resistance, fits well to the case of pure but not totally defect-free metal 



samples in the whole temperature range measured [40], as shown in Fig. 3.9 The absolute 
value of the Wiedemenn-Franz ratio is found to be ~ 2.53 x W~ s W£l/K 2 at 290 K, slightly 
larger than that of the free-electron Lorenz number 2.45 x 10~ 8 W$7/K 2 , which is also in 
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Figure 3.8: Experimental test of the 3o> method on Pt wire of 20 /jm in diameter, (a) 
Current dependence of the V^. The open circles are the measured data, and the solid line 
is the predicted relation oc I 3 , (b) Frequency dependence of the (open circles). 
The solid line is the predicted relation oc 1/ ' ^Jl + (2co^) 2 . (c) Frequency dependence 
of the phase angle of the (open circles). The solid line is again the predicted relation 
tanc/> oc wy. The thermal conductivity k, thermal diffusivity a 2 , and specifc heat C p of the 
Pt specimen are plotted as open circles in figures (d), (e), and (f), respectively. Also shown 
as solid circles in these figures are the standard data of Pt [39]. The difference in k and a 2 
between our measured data and the standard data should reflect the difference in purity 
and/or structural perfection between the Pt specimens of different sources. 
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T/e 

Figure 3.9: The Wiedemenn-Franz ratio L e of our Pt specimen compared with that of usual 
metals with different purity. The result indicates that our Pt wire specimen is pure but not 
totally defect-free. The room-temperature Wiedemenn-Franz ratio of our Pt specimen is 
about 2.52 x 10 _8 WS7/K 2 , in perfect agreement with what reported in literature [40], but is 
slightly larger than that of the free-electron Lorenz number Lq = 2.45 x 10 _8 Wil/K 2 (the 
dashed line). Note that for Pt its Debye temperature 6 is 240 K. 



perfect agreement with the previously reported value in literature [40] . 

Let us now examine the effect of radial heat loss through air conduction. Different 



sets of data in Figs. 3.10 (a), (b) and (c) were taken at different vacuum pressures. Open 
squares correspond to the data taking in high vacuum where slight changes (e.g., a factor 
of 2) in vacuum pressure does not affect the results, i.e., heat loss through air is absent. 
In poor vacuum, however, such heat loss results in spurious larger thermal conductivity 



and diffusivity of the specimen (open circles in Figs. 3.10 (a) and (b)). Nevertheless, it is 
interesting that the specific heat deduced from the k and a 2 data is relatively insensitive to 



such heat loss (Figs. 3.10 (c)). 



After all, let us check if the approximation conditions (3.13) and (3.33) are satisfied. 
llR'L 



We find that gr~jr~g ~ 10 (taking n = 1). Therefore, (3.13) is well satisfied. For radiation 
heat loss, we find that g ~ 0.22 s _1 even at 300 K, and 7 (actually, 7') deduced from the 
measurement is about 0.2 s. Therefore, 57 ~ 0.044, which means that omitting the radiation 
heat loss will cause about 5% of error for a 20 ^m diameter Pt wire specimen with length 
of 8 mm. 

We have also applied the 3a; method to measure the k and C p of multiwall carbon 
nanotube bundles who have a negative temperature coefficient in electrical resistivity. The 
bundles are about ~1 mm long and with a diameter ranging from ~1 /xm up to ~100 
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Figure 3.10: Effect of radial heat loss through air conduction. Open squares represent the 
case where such heat loss is absent. In one run of the measurement, radial heat loss was 
triggered on above T* (marked in the figure) by distroying the vaccum of the system. The 
heat loss results in spurious larger thermal conductivity and diffusivity of the specimen 
(open circles in (a) and (b)). Nevertheless, the specific heat deduced from the k and a 2 
data is relatively insensitive to such heat loss (Fig. (c)). 
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/im. For these specimens, no C p and k data of other sources are available for comparison. 
Nevertheless, the frequency dependencies of V% u and tan0 obtained on the specimens are 



all in good agreement with the predicted functional forms (see Fig. 3.1), which guarantees 



the reliability of the deduced k and C p . For the specimens in this measurement, ~^r^ is 
below 10~ 3 above ~ 60 K and is about 0.08 at 10 K, and <?7 is below 5xl0~ 3 in the whole 



temperature range, both of which satisfy the conditions (3.13) and (3.33), respectively. For 
a 1 fim diametered specimen, the mass is about 10 -9 g, which is far less than the minimum 
amount of mass required in any other kinds of C p measurement. 

3.5.3 Summary 

We have presented an improved 3w method after explicitly solved the heat-conduction 
equation under some reasonable approximations. This method can simultaneously measure 
the thermal conductivity and the specific heat of electrically conducting materials in a way 
similar to, and as simple as usual electrical resistivity measurement. 
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Chapter 4 

Tunneling Spectroscopy 



4.1 Motivation 

Coulomb blockade (CB) has been intensively studied in a multij unction configura- 
tion, in which electron tunnel rates from the environment to a capacitively isolated "island" 
are blocked by the e-e interaction if the thermal fluctuation is below the charging energy 
E c = e 2 /2C and the quantum fluctuation is suppressed with sufficiently large tunnel resis- 
tance Rt 3> Rq = h/2e 2 . In the case of a single-junction circuit, the understanding of CB is 
less straightforward. The Coulomb gap, supposed washed out for the case of low-impedance 
environments, should only be established if the environmental impendence exceeds Rq. 

In singlewalled carbon nanotubes (SWNTs), the reduced geometry gives rise to 
strong e-e interaction. Indeed, CB oscillations and evidence of Luttinger liquid (LL) have 
been observed [1]. In contrast to SWNTs, in which only two conductance channels are 
available for current transport, MWNTs with diameter in the range of 20-40 nm have 
several tens of conductance channels. The energy separation of the quantized subbands, 
given by AE = hvf/d, is about 13-26 meV taken Fermi velocity as Vf = 8 x 10 5 m/s. 
This value is about an order of magnitude smaller than that of SWNTs. Experiments 
indicate that MWNTs are considerably hole doped, thus a large number of subbands in 
the order of ten are occupied. The observation of weak localization [2], electron phase 
interference effect [3] and universal conductance fluctuations [2] support the view that low 
frequency conductance in MWNTs is contributed mostly by the outmost graphene shell 
and is characterized by 2D diffusive transport. Besides the phase interference correlated 
effects, strong e-e interaction has also been observed in MWNTs. There have been several 
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tunneling spectroscopy observations of a pronounced zero-bias suppression of tunneling 
density of states (TDOS) [4-6]. Moreover, the TDOS shows a power law asymptotics, i.e., 
v{E) ~ E a , which resembles the case of a LL. It is noteworthy that in the Environmental 
Quantum Fluctuation (EQF) theory, for a single tunnel junction coupled to high-impedance 
transmission lines, such a scaling behavior is also predicted in the limit of many parallel 
transmission modes [7]. The physical origin of these power laws is the linear dispersion 
of bosonic excitations that are characteristic both for LL, which is strictly ID ballistic 
conductor, and a single tunnel junction connected to a 3D disordered conductor. In the 
latter case the quasi particle tunneling is suppressed at V <C e/2C, therefore the charge 
is transported with ID plasmon modes. The facts that MWNTs have many conductance 
modes, as well as the observation of a crossover from power law to Ohmic behavior at higher 
voltages [5], suggests that EQF theory is more appropriate to describe the observed TDOS 
renormalization. Most of these measurements were done in multi-junction configurations. 
Therefore a single tunnel junction measurement is needed to further clarify this issue. 

4.2 Experimental 

In this chapter, millimeter-long CVD-grown MWNTs [8] are measured by a cross- 
junction method. The MWNTs samples are composed of loosely entangled nanotubes that 
are roughly parallel to each other and extend up to 2 mm long. The single tunnel junctions 
are formed by crossing a very thin (< 1 jum) MWNTs bundle on a narrow strip of metal wire 
fabricated on an insulating substrate. We explored different electrode materials including 
Au, Cu, Sn, and Al. In the case of Sn and Al, a small magnetic field was applied to suppress 
the Superconducting state below Tc. No obvious change of device characteristics caused by 
the choice of metals was observed. With such a cross-junction configuration the measured 
conductance is only contributed by the tunnel junction itself, since the current is passed 
along one arm of the MWNT/metal and voltage is measured along the other non-current 
carrying arm. Thus the device can be understood as a small number of single junctions in 
parallel. Despite the simplicity of the fabrication method, we found that the devices made 
are very stable and sustain several cooling cycles without apparent change of characteristics. 
More than 20 samples we measured all yield strong zero bias suppression of TDOS. 

The current I(V) as well as G = dl/dV, are calculated by a golden-rule approach 
incorporating the environmental influence by P{E), the probability for a tunneling electron 
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V/(e/2C) K > 

Figure 4.1: (a) dl/dV as a function of V in dimensionless units measured at T = 0.12, 0.25, 
0.78, 1.2, 1.7, 4.5, 10, 20 K (sample 990530s6). Dashed lines are the fits of EQF theory, (b) 
DC current simultaneously measured at 0.12 K. Dashed line is the numerical fit. Dotted 
line shows a Coulomb offset of e/2C. (c) Temperature dependence of the exponent a for 
three samples (O 990530s6; A: 990316; y: 990320). 
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to lose energy to the environment. P(E) can be calculated from its Fourier transform, the 
phase correlation function J(t). In the transmission- line model, J(t) is determined by the 
total environmental impedence 

Z t (u) = (iuC + Z- 1 (u))- 1 (4.1) 

where Z{oS) is the external environmental impedance. We applied the method in Ref. [9] 
to evaluate P{E) from an integral equation without the need of going to the time domain. 
The I-V characteristics are then calculated for finite temperatures and arbitrary junction 
impedance. The parameters that appear in the numerical calculation are the damping 
strength 

a = Z(0)/R Q , (4.2) 
the inverse relaxation time wrc, an d the quality factor 

Q = u R c/us (4.3) 

with the resonance frequency of the undamped circuit us = (LC)^ 1 ^ 2 . It is noted that the 
quality factor Q only plays a minor role as an additional adjustable parameter. Therefore 
it is not important to the fit, and is set always to unit. 



4.3 Discussion 

We then tried to fit our experimental data with the EQF theory. Since E c and Rt 
are determined by the high-voltage data, the isolation resistance Ri so = Z(Q) is the only 
adjustable parameter in our fitting. In contrast to the case of Ref. [10], where Ri so is formed 
by ideal Ohmic and temperature-independent resistors, Ri so in our case are provided by the 
resistive impedance of the MWNTs themselves, which should be temperature-dependent. 
From the fitting, we do see an evolution of a from 0.1 at 20 K to 0.25 at 1-4 K and then 
a saturates (Fig. 4.1;). Note that dimensionless units are used with G normalized by 1/Rt 



and voltage normalized by e/2C. Therefore the number of the single tunnel junctions has 
no effect on the fitting and results on different samples can be directly compared. We found 



that for different samples with scattered characteristics (see Table 4.1 ), at low temperatures 
the exponent a reaches an universal value of 0.25-0.35, which agrees with previous multi- 
junction measurements [5,6]. It yields Ri so = 3.3-4.6 kO that is roughly a constant for 
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Table 4.1: A partial list of characteristics of the samples. 



Sample 


Rt (kfi) 


E c (eV) 


C (aF) 


Rise (kfi) 


a 


990316 


18.4 


0.019 


4.2 


3.3 a 


0.25 


990320 


48.3 


0.014 


5.7 


3.3 a 


0.25 


990530s6 


8.6 


0.01 


8.0 


3.3 a 


0.25 


990530sl 


10.0 


0.02 


4.0 


4.6 a 


0.35 


990202 


9.1 


0.004 


20.0 


3.5 a 


0.27 



a Value acquired at T = 1.2-4.3 K. 



different samples. This coincidence is not accidental but reflects the intrinsic electro dynamic 
modes of the MWNTs, which can be modeled as an ideal resistive LC transmission line: 

R = (L'/C) 1 / 2 (4.4) 

with the kinetic inductance L' = Rq/2Nvf estimated as ~ 1 nH/^m for iV ~ 10-20 modes 
and the capacitance C ~ 20-30 aF//um. Therefore, the "environment" with respect to the 
single-junctions in our devices is provided by MWNTs themselves, not the external circuits. 

As mentioned previously, for an isolated single tunnel junction with many parallel 
transport modes, the EQF theory predicts a power-law asymptotics: 

G «T)/ G (0,T) S HV/T) = |r[f + ^]/r[f Ml + ^]|> (4.5) 

where G(0,T) is the zero-bias conductance [10]. For eV/2irkBT 3> 1, a voltage power law 
G(V, T) ~ V a is expected. Note that in the above scaling function, exponent a is the 
only adjustable parameter. In a two-junction configuration, an additional parameter was 
introduced to take into account the voltage division between the two tunnel junctions [1]. 



Seen in Fig. 4.2, we do observe such a scaling behavior. Here f(V/T) represented by solid 
line is calculated taken a = 0.27. The inset shows the original data in which the dashed 
line is calculated numerically using the same exponent a = 0.27. 

CB and Coulomb gap effects in MWNTs are recently described with microscopic 
theories in Refs. [14,15]. In the 2D diffusive regime, the exponent for tunneling into the 
bulk of a MWNT is 

a = (R/hu D)ln(l + u U ) (4.6) 

Here R is the tube radius, Uq is the intratube Coulomb interaction, D = Vjt/2 is the charge 
diffusivity, and the "bare" DOS vq = N/Ahvf with N ~ 20. If we use vtr ~ 60 nm from the 



Chapter 4: Tunneling Spectroscopy 



58 




Figure 4.2: Scaled conductance G(V,T)/G(0,T) of another sample (990202). The inset 
shows the original dl/dV vs V with the dashed line calculated with exponent a = 0.27. 



Chapter 4-' Tunneling Spectroscopy 



59 



magnetoresistance data [11] and taking R ~ 20 nm, then under the condition of Uo/hvf ~ 1 
the above formula yields a ~ 0.25, which agrees with our experiment. 

4.4 Fano resonance 

Besides the voltage power-law phenomenon, which is characteristic of strong Coulomb 
interactions, another energy scale "the discrete energy levels due to electron confinement" 
emerges as temperature decreases, and it modifies the tunneling spectra. Unlike the case 
of a two-junction configuration, where the electrons are confined by the two contacts if the 
nanotube is clean enough, in a single-junction configuration the electrons can be confined by 
disorders when the MWNTs are "dirty" , such that the impedance of the local environment 
of the junction is larger than Rq. A stacking mismatch between adjacent walls and other 
structural imperfections are possible sources of disorders in MWNTs, resulting in discrete 
energy levels. 

By cooling down the devices to even lower temperatures, we observe that G devel- 



ops a narrow resonance- like anomaly at very low bias (Fig. 4.3). The asymmetric anomaly 
builds up consistently with the decrease of temperature and even shows a dip structure. 
The line shape resembles that of a Fano resonance, which is ubiquitous in resonant spin-flip 
scattering. Fano resonance happens when there are two interfering scattering channels: a 
discrete energy level and a continuum band. It has recently been rediscovered in mesoscopic 
systems such as semiconductor quantum dots, SWNTs, etc. [12-14]. Taking a quantum dot 
as an example, it can be considered as a gate-confined droplet of electrons with local- 
ized states. The coupling of the dot to the leads can be tuned in a controllable way to 
enter different transport regimes: If the dot is weakly coupled from the environment, a 
well-established CB is developed. The charge transport is suppressed except for narrow res- 
onances at charge degeneracy points. When the tunnel barriers become more transparent, 
the dot enters Kondo regime, i.e., below a characteristic Kondo temperature Tk, spin-flip 
co-tunneling events introduce a narrow symmetric TDOS peak at Ep that can be inter- 
preted as a discrete level. If the coupling is strong enough, the interference between this 
discrete level and the conduction continuum give rise to an asymmetric Fano resonance. 
In Ref. [12], such a crossover from a well-established CB through Kondo regime to Fano 
regime has been clearly observed. 

We found that the asymmetric resonance line shape of G can be fitted by the 
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Figure 4.3: Asymmetric resonance features of dl/dV seen in two samples: 990530sl(a) and 
990530s6(b). Dashed lines are the fits of Fano's formula. 



Table 4.2: Characteristics of the two samples in Fig. 4.3 
Sample T (K) q e (meV) T (meV) T K (K) 



990530s6 


0.25 


0.53 


0.43 


0.32 


1.86 


990530s6 


0.12 


1.75 


0.57 


0.29 


1.68 


990530sl 


0.25 


0.92 


0.3 


0.65 


3.77 


990530s! 


0.08 


1.97 


0.43 


0.56 


3.25 



Fano's formula: 

G oc (E + qf/{E 2 + 1) (4.7) 

Here q is the so-called asymmetry parameter. E = (e — eo)/(F/2) is the dimensionless de- 
tuning from resonance. As expected, the asymmetric parameter q as a measure of the degree 
of coupling between the discrete state and the continuum increases when the temperature 
drops. The Kondo temperature, estimated from the relation V = 2ksTk, is consistent for a 
device at different temperatures, and agrees with the observed FWHM of the resonance. 

If Kondo physics really accounts, then G at voltage near resonance should show a 
non-monotonic temperature dependence around Tx [14]. A close look at I-V curves (Fig. 



4.3) does exhibit such a behavior: G at the peak position first drops with T and then rise 



up below ~ 1 K. The conductance measured at zero bias, however, shows a monotonic 
temperature dependence (Fig. |4.4b ), which is expected when measuring at energies out 
of the Fano resonance. Moreover, we found that a perpendicular magnetic field B shows 



two- folded effects (Fig. 4.4 1 and 4.4;): Firstly the background conductance monotonically 
increases with the applied B. Secondly, the dip seen at zero field gradually disappears and 
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T (K) V (mV) 

Figure 4.4: (a) The dependence of zero-bias dl/dV on perpendicular magnetic field, (b) 
The temperature dependence of dl/dV measured at different voltages, (c) Effect of per- 
pendicular magnetic field. From bottom to top: B = 0, 1,4, 8 Tesla (sample 990530s6). 

the resonance turns into a nearly symmetric peak at high field, similar to what was observed 
in the semiconductor quantum dots [12]. The magnetic field effect can be explained with 
the following pictures: Firstly adding a flux should change the amplitude and/or the phase 
for the resonant channels therefore breakdown the coherent backscattering and increase the 
forward transmission through the channel. Secondly, the magnetic field can destroy the 
interference between the resonant and nonresonant path, transforming a resonant dip into 
a peak. 

Since Kondo physics is historically intepreted as the interplay between the d- 
orbitals of magnetic impurities and the conduction continuum, we have to exclude the 
possibility that the observed Fano resonance come from the residual Fe/Si catalyst in the 
MWNTs samples. It has been shown that magnetic impurities in MWNTs cause an enhance- 
ment of thermoelectric power [15], which is absent in our careful TEP measurement [16]. 
Furthermore, no Fe signature in the body of MWNT bundles could be detected within the 
instrumental resolution in the X-ray and TEM studies. Therefore the observed Fano reso- 
nance should be inherent property of MWNTs, which can be virtually treated as quantum 
dots strongly coupled to the leads. We noticed that recently Buitelaar et al. have also 
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seen clear traces of CB and Kondo resonance in MWNT SET devices [17]. Their observed 
Kondo temperature Tk = 1.2 K is in good agreement of our results. 

4.5 Summary 

In summary, strong zero-bias suppression of TDOS is observed in MWNTs single 
tunnel junctions that can be well explained with the EQF theory. The observed exponent 
a ps 0.25-0.35 is found to be consistent for all the samples. Reinforced by the observation of 
Altshuler-Aronov-Spivak (AAS) oscillations in magnetoresistance measurements (presented 
in the next chapter), we interpret MWNTs as a quasi-2D diffusive molecular quantum wire 
in which both disorder and intrinsic Coulomb interaction play important role in low-energy 
excitations. Similar to the case of an open quantum dot, in low-impedance junctions we 
found that a Fano-resonance like asymmetric conductance anomaly is built up below meV 
energy scales. The coexistence of abundant mesoscopic phenomena exhibits that MWNTs 
provide a good laboratory for the study of electron interactions in mescoscopic system. 
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Chapter 5 

Transport Properties 



5.1 Thermoelectric power 

5.1.1 Introduction 

The thermoelectric power is an important transport property of solids. The basic 
picture of thermoelectric effect (or Seebeck effect) in a conductive solid is the carrier diffusion 
in the presence of a temperature gradient. The specimen is considered as a vessel filled with 
a gas of free electrons and/or holes, and when a temperature gradient A X T = dT/dx is 
applied, conduction electrons and/or holes are thermally drifted in the -x direction, so as 
to setup an electric field E x as well as a thermoelectromotive force (TEF) 

V x = I E x dx (5.1) 
J o 

The thermoelectric power (TEP) or Seebeck coefficient is defined by 

In the case of a two-carrier system, electrons and holes give rise to polarizations 
opposite to each other, so that S < when electrons predominate over holes in the contri- 
bution, and vice versa. 

In the above discussion, the internal structure characters of materials are not taken 
into account. In actual solids, thermal vibrations of atoms can be described as combined 
elastic lattice wave, or in the quantum mechanical expression, these lattice waves are con- 
sidered as phonons. In absense of electron-phonon interactions, a specimen can be looked 
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on as a vessel of electron and phonon gases in which two constituents diffuse separately in 
the presence of a temperature gradient. However, in actual solids the momentum of the 
drifting carriers is irreversibly transfered to phonons and give rise to electrical resistance 
and Joule heating. In the thermoelectric effect, since the drift takes place in the same direc- 
tion for both systems, phonons may transfer their momentum to carriers so as to enhance 
the polarization. Such a mechanism, the so-called phonon drag effect, has been manifested 
in a variety of solid materials. 

Since carbon nanotube can be considered as wrapped-up graphite sheets, it is 
helpful to discuss the thermoelectric power of graphite firstly. The well crystallized graphite 
is a nearly compensated material, therefore the diffusive component of TEP is very low and 
its sign changes very sensitively, reflecting the delicate balance between electrons and holes. 
When the electronic structure is concerned, graphite is a semimetal in which the phonon 
drag effect is pronouced due to the following reasons. Firstly, the electron-phonon coupling 
which participates directly in the drag effect is very strong. Secondly, the wave vectors 
of the phonons, which play a dominant role in the drag effect, are small due to the small 
Fermi surface of graphite. Therefore the phonon-phonon collisions establishing equilibrium 
phonon distribution is not strong enough to diminish the phonon drag effect. As a result, 
the basal-plane TEP of graphite has a pronounced phonon drag induced dip at about 30K, 
at which the average wavenumber of the in-plane mode phonons is found to equal the Fermi 
wavenumber. 



5.1.2 Experimental 

We use a standard comparative method to measure the TEP and thermal con- 
ductivity of a MWNT bundle sample simultaneously. The measurement configuration is 



shown in Fig. 5.1. The sample was placed in series with a standard constantan (TN) rod 
between the electric heater and the heat sink. Two chromel (KP) wires were attached to the 
constantan rod to form a thermocouple. Assuming that there were no heat losses, i.e., the 
same heat flow passes through the constantan rod and the sample. Four Au electrical leads 
were attached to the sample with silver paint so that successive measurement of electrical 
resistance could be performed on the same sample. Chromel-constantan thermocouple leads 
were attached with high thermal-conductive epoxy at the same place with the voltage leads. 
TEP and thermal conductance was measured at fixed temperatures by heating the far end 
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Heater 




Figure 5.1: Schematic of experimental configuration for TEP measurement. The two Au 
current leads for four-probe resistance measurement are not shown. 

of the constantan rod with triangle- wave formed currents at very low frequencies (i.e., 0.01 
Hz). By comparing different voltages induced on the constantan rod and the sample, we 
obtained the thermal conductance of the tubes. 

According to the formula of Seebeck effect, S of two chromel-constantan thermo- 
couples and MWNT sample are as follows 

V x = (S TN - S K p)AT t (5.3) 
V 2 = (S TN - S KP )AT 2 (5.4) 
V3 = {Smwnt - Sau)AT 2 (5.5) 
The heat currents flow in the constantan rod and MWNT bundle could be written 

as 



Qi = -ktn-AiATi/Li 
Q2 = —KMWNTA2AT2/L2 



(5.6) 
(5.7) 
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Figure 5.2: TEP of a MWNTs sample measured from 10 K to 300 K. 

Here A\ and A2 are cross-sectional area of the constantan rod and MWNT bundle 
separately. L\ and L2 are their length. Assuming that there were no heat losses at the 
interface between the constantan rod and the MWNT sample, i.e., Q\ = Q2, S and k of 
MWNT sample can be derived 



Smwnt = Sau + (Stn — Skp)Vz/V2 
kmwnt = ktn{AiL 2 / A 2 Li)Vi/V2 



(5.8) 
(5.9) 



Shown in Fig. 5.2 are the results of our measurement. TEP of MWNTs is mod- 



erately large, approaching ~ 20 //V/K at room temperature. The sign of TEP is positive, 
indicating that holes predominate over electrons in the contribution. At low temperatures, 
TEP has a nearly linear temperature dependence and approaches zero as T — > 0, which is a 
character of metallic behaviors. At high temperatures TEP tends to be saturated. Such a 
large positive thermoelectric power of MWNTs gives strong evidence that the electron-hole 
symmetry in metallic nanotubes is broken. Similar results on SWNT ropes were reported 
by Hone et al. [1]. At the present stage, the origin of such symmetry breaking is unclear. 
Hone et al. suggested that the electron-hole symmetry is broken due to inter-tube hopping 
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when the tubes are assembled into ropes. However, our MWNTs are well decoupled from 
each other, therefore such a mechanism is not appropriate. Other mechanisms, such as 
carrier interactions with disorder and defects, might play important roles in TEP. 

The data of thermal conductivity (not presented) is qualitatively similar to that 
derived by 3w method. In most of the temperature range measured, n is nearly linear. The 
order of magnitude of k is the same with that derived by 3a; method, giving a convincing 
proof that 3w method is very reliable. 



5.2 Magnetoresistance 

In this section, the dependence of the electrical resistance of MWNTs on longitu- 
dinal magnetic field (parallel to tube axes) has been studied. Without access of lithographic 
micro-fabrication facilities, we prepared the devices by the following method: A thin bundle 
of MWNTs was laid on a sapphire substrate with parallel chromium electrodes separated 
with each other by 4 fj,m. Then a zircite wafer with a layer of single-sided Scotch tape 
was pressed on the sample by a tiny screw. The Scotch tape was used to prevent the 
breaking of the sample. By adjusting the magnitude of pressure, good electrical contacts 
could often be made. The room-temperature two-probe resistance of the samples range 
from several tens kf2 to 1 MO. The experiment was conducted in a dilution refrigerator 



with a superconducting magnet (Oxford Instruments 200 TLM). Figure 5.3 shows a typical 
longitudinal magnetoresistance (MR) measurement. When a magnetic field B is applied, 
the resistance of MWNTs decreases. This negative MR is associated with the phenomenon 
of weak localization (WL). WL originates from the enhancement of backscattering which 
contains interference terms adding up constructively in zero field. In a specific geometry 
like a diffusive and thin- walled metallic cylinder (or ring), the WL contribution is periodic 
in the magnetic flux through the cylinder with period h/2e, known as Altshuler-Aronov- 



Spivak (AAS) effect [2,3]. In Fig. 5.3 it can be clearly seen that the resistance has a second 
maximum at B = 3 T and a third maximum at B = 6 T. From TEM images, the average 
outermost diameter of our MWNTs is 30 nm, the corresponding period of AAS oscillations 
is 2.93 T, which is in good agreement with the observation. 

Besides the reproducible resistance oscillations, non-reproducible temporal resis- 



tance jumps are often observed (Fig. 5.4). We note that Bachtold et al. observed similar 



jumps in measurements on individual MWNTs. They suggested that temporal changes in 
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Figure 5.3: Longitudinal Magnetoresistance of a MWNTs sample at 20 mK (solid curve). 
The temporal resistance jumps of the bare data have been removed. The dashed curve 
corresponds to theoretical prediction for the quantum correction to the resistance of a thin- 
walled cylindrical conductor using the parameters mentioned in the text. 

the electric contacts, defects or impurities in the nanotubes are possibly responsible for 
these jumps. However, we note that the magnitudes of these temporal resistance jumps 
are integer multiples of h/2e 2 . Therefore, the resistance jumps probably correspond to 
the openning or closing of the conducting channels of a MWNT. Further investigations are 
needed to clarify this problem. 

Now let's compare our experimental data with the theory of AAS effect [2,4]. 
Firstly we consider a single graphene cylinder with radius R and longitudinal height 6, 
which approximately equals to the distance between the electrical contacts. Then two 
parameters appear in the AAS theory: the temperature-dependent phase-coherence length 
Lip and the mean free path I. Following the formula given by Altshuler, Aronov, and Spivak, 
the correction of longitudinal DC conductance is 

_ e 2 2nR„ Lin ^ . 27i\R N , 26 x , 

AG = 5--— In -f + 2 V K (n—) cos 27rn-f 5.10 

Tr z h b I ^ Lip <po 

n=l r 

where K${x) is the modified Bessel function. The conductance correction oscillates 
with a period <fio/2, becoming exponentially small for L^ <C 2ttR. 
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Figure 5.4: The raw data of the same sample as in Fig. 5.3 Temporal resistance jumps 
with magnitude h/2e 2 , h/e 2 and 2h/e 2 can be clearly seen. 



In Eq. |5. 10 the wall thickness was not taken into account. In the real samples 
the field in the cylinder walls is not negligible. As a result, the oscillations decay as the 
field increases and a monotonic component appears in the magnetic field dependence of 
the sample resistivity. Moreover, the magnetic field direction may deviate from the axis of 
cylinder by a small angle 9, making the effective wall thickness a* larger than the actual 
value a. The inclusion of magnetic field results in the phase relaxation length L„ becoming 
field dependent. For a cylinder with an effective wall thickness a* in a parallel field, 

1 1 l,a*efl\. 



L 2 Dt w + 3^ n 



where a* is given by 



. *2 



a 2 cos 6 2 + 6R 2 sin 9 2 



(5.11) 



(5.12) 



Matching the theory to the measurement at T = 20 mK, we obtain L v ~ 63 nm, / ~ 1 nm, 
6 = 4.4°. The agreement between theory (dotted line) and experiment (solid line) is quite 
good except for the low field part. This agreement is obtained when assuming that only one 
graphene cylinder contributes to the conductance, i.e., the wall thickness a = 0.34 nm. The 
estimated value of the phase relaxation length L„ is comparable with the circumference of 
graphene cylinder, i.e., L v « 2nR, and is close to that derived by Bachtold et al. [5] on 
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individual MWNT at T = 0.3 K. Theory also points out that the phase of the oscillations, 
i.e., m&xR for H = 0, will be reversed by the presence of spin-orbit scattering. In our case, 
phase reversal does not happen, indicating that spin-orbit scattering is negligible. 

Therefore, we conclude that quantum-interference corrections to the resistance 
can account for the measured longitudinal MR. The good agreement between theory and 
experiment give a convincing proof that the electric current flows in the outermost metallic 
graphene wall of MWNTs at low temperatures. 

5.3 Summary 

Transport properties, including thermoelectric power (TEP) and longitudinal mag- 
netoresistance (MR), have been investigated on millimeter- long aligned multiwall carbon 
nanotubes (MWNTs). A large and positive TEP with a nearly linear temperature de- 
pendence is observed, indicating that the electron-hole symmetry in metallic nanotubes is 
broken. The low-temperature longitudinal MR of MWNTs shows an oscillation behavior, 
which is attributed to Altshuler-Aronov-Spivak (AAS) effect. This gives clear evidence 
that quantum interference phenomena dominate the MR at low temperatures. The phase 
relaxation length is found to be comparable with the circumference of the outmost 
graphene wall. The observation of resistance jumps at integer multiples of h/2e 2 suggests 
that MWNTs are mesoscopic conductors with finite conduction channels. 
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Conclusions 

Main conclusions from our studies on the thermal and electrical properties of 
MWNTs are as follows: 

1. The phonon structure of a MWNT with diameter of several tens of nm is two- 
dimensional. The inter-wall interactions in a MWNT is very weak, so that one can treat a 
MWNT as a few decoupled two-dimensional (2D) single wall tubules. 

2. The low thermal conductivity indicates that large amounts of defects exist in 
the MWNTs derived by thermal decomposition method. 

3. Strong Coulomb blockade induced zero-bias suppression of tunneling density 
of states (TDOS) is observed in single tunnel junctions between MWNTs and a normal 
metal. Quantitative analysis by Environmental Quantum Fluctuation (EQF) theory found 
excellent agreement with experimental data. A universal exponent a ~ 0.25-0.35 is at- 
tributed to strong intratube Coulomb interaction in MWNTs. In low-impedance junctions, 
Fano-resonance like asymmetric conductance anomaly is found below mV energy scales. 

4. The thermoelectric power (TEP) of MWNTs is moderately large. At low 
temperatures, TEP has a metallic-like linear temperature dependence. Such a large positive 
thermoelectric power of MWNTs gives strong evidence that the electron-hole symmetry in 
metallic nanotubes is broken. 

5. The longitudinal magnetoresistance (MR) of MWNTs shows an Altshuler- 
Aronov-Spivak (AAS) effect induced oscillation at 20 mK. This result gives clear evidence 
that quantum-interference dominate the resistance of MWNTs at low temperatures. The 
estimated value of the phase relaxation length L v is comparable with the circumference of 
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the graphene cylinder. 
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This dissertation is revised from the original version submitted seven years ago. 
At the time I had not learned how to include figures in LaTeX. Instead, figures were hard 
copied and the thesis was submitted in printed format only. The main purpose of the 
revision is to reserve a digital copy and make it available online for interested readers. For 
authenticity, the main content of the thesis is kept As Is, except for the necessary corrections 
and reformatting. Noticeable changes are listed as follows: 

1. Chapter 4 (Tunneling Spectroscopy) is updated according to the recent publication 
(PRL 91, 076801 (2003)). The corresponding part in Abstract is updated as well. 

2. In Chapter 1, Fig. 1.3(b) is changed. Original figure (curves of constant energy of 
graphene sheet) is changed to the current version (DOS of carbon nanotubes from 
STM and tight-binding calculations). Fig. 1.5 is corrected. Original figure (band 
structure of (5, 5) nanotubes) is replaced by band structure of (10, 10) nanotubes to 
make it consistent with the context. 

3. In Chapter 2, two new figures (Fig. 2.1 and 2.4) are added as complementary mate- 
rials. Fig. 2.3 is replaced by the published version from Ref. [3]. 

4. In Chapter 3, Fig. 3.3 (from Ref. [22]) is replaced by more recent data (from Ref. 
[17]). Fig. 3.5 and 3.9 (omitted in the original thesis) are now included. 

5. In Chapter 5, Fig. 5.1 (configuration for TEP measurement) is redrawn to include 
more details. 

6. Summary sections are added for Chapter 2 and Chapter 5. 

7. English translation of Acknowledgments is added. 

A recent survey by Michael Banks at the Max Planck Institute has ranked carbon 
nanotubes as one of the top five topics in physics publications. Although this does not 
promise that someday they would replace silicon microelectronics, or make space elevators, 
carbon nanotubes and related materials remain attractive for both fundamental research 
and potential applications. Since the nanotube research is evolving very fast, the content 
of this thesis may become obsolete. Somehow I still feel obligated to finish the revision, not 
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just for fun, but also to make a repository of the pleasant memories of my first scientific 
safari. 
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